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£^ . We investigate the scaling properties of the specific heat of the XY model 

7— I \ 

on lattices H x H x L with L 3> H (i.e. in a bar-like geometry) with respect 

q . to the thickness H of the bar, using the Cluster Monte Carlo method. We 

<D ; 

study the effect of the geometry and boundary conditions on the shape of 

+-> , 

the universal scaling function of the specific heat by comparing the scaling 

(Z3 ; 

-*-J . functions obtained for cubic, film, and bar-like geometry. In the presence of 

physical boundary conditions applied along the sides of the bars we find good 

'■ 

agreement between our Monte Carlo results and the most recent experimental 

o : 

data for superfluid helium confined in pores. 
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' Bulk liquid 4 He exhibits a second order phase transition at the A-critical temperature 

<n : 

T\ ~ 2.18K where it turns superfluid. Approaching T\ from below the specific heat c and 

OO ' 

the superfluid density have singularities which are characterized by the critical exponents 
a and u, respectively; these critical indices determine the universality class of this second 
order phase transition. If liquid 4 He is placed in a confining geometry, e.g. film or pore- 
like geometry, the finite-size scaling theory || can be used to describe the behavior of the 
physical quantities at temperatures close to T\. The finite-size scaling theory is based on 

•i-H , 

the assumption that the system feels its finite size when the correlation length £ becomes of 
the order of the confining length. For a physical quantity O this statement can be expressed 
as follows 0: 

0(t,H = oo) J[h [) 

1 



where H denotes the relevant confining length and the reduced temperature t = T/T\ — 1 
while £(t,H = oo) is the correlation length of the bulk system. The point is that the 
dimensionless function / depends only on the dimensionless ratio H/£ and it does not 
depend on microscopic details of the system. It does, however, depend on the observable O, 
the type of confining geometry and on the conditions imposed (or not, in the case of free 
boundaries) at the boundaries of the system. 

Liquid 4 He can be an ideal testing ground to check the validity of the finite-size scaling 
theory experimentally || because the specific heat c and the superfluid density can be 
measured to a very high accuracy. In addition, the shape of the confining geometry, such 
as films or pores, can be designed with such a precision that the relevant confining length is 
well defined. 

In order to compare the results of finite-size scaling theory with the experimental results 
on liquid 4 He, there have been theoretical efforts to compute the universal finite-size scaling 
function of the specific heat of confined liquid 4 He PHTOf. The scaling function for film 
geometry has been computed by loop expansion based renormalization group treatment of 
the standard 4 Landau- Ginzburg functional an d by the Monte-Carlo method within 



the XY model P|-|TT](]. The XY model is another form of the standard Landau- Ginzburg 
functional and both models belong to the same universality class as liquid 4 He, i.e. all three 
systems have the same critical indices. The scaling functions for a film geometry have been 



computed from the XY model using the Monte Carlo method |q-[T0| and they are in rather 
good agreement with the scaling functions obtained from the field theoretical treatment of 
the 4 Landau- Ginzburg theory ||. The scaling functions obtained when periodic boundary 
conditions |J were applied on the top and on the bottom boundary of the film are very 



different from the scaling functions obtained when staggered boundary conditions p|JT0 
(vanishing order parameter on the boundary layer) were applied on the top and bottom film 
boundary. The scaling functions obtained by the Monte Carlo simulation of the XY model 
PJiTl|] agree reasonably well with the experimentally determined || scaling functions when 
staggered boundary conditions were applied in the simulation. Such boundary conditions 



on the order parameter are believed to approximate better the real conditions on the top 
and bottom boundary of superfluid film |6].|8| . 

Besides the boundary conditions, the other important factor which, in principle, can 
determine the scaling function is the geometry. Close enough to the critical point where 
the bulk correlation length £ becomes of the size of the confining length H or even larger, 
the system as a whole "knows" about the geometrical shape of the confining volume. Thus, 
one expects to find a region of the dimensionless variable H/ £ where the form of the scaling 
function is sensitive to the geometry of the confining volume. There are experiments which 



have been performed for the pore geometry [|TTH T3|1 and the results for the specific scaling 



function are significantly different from those obtained from the film geometry. In this 
paper, we report results of our simulations of the XY model in a bar-like geometry, namely 
on lattices H x H x L with L » H, which mimic the pore-like geometry of the experiments. 
In addition, to represent the real situation more closely we have used staggered boundary 
conditions along the confining dimensions (H directions). Periodic boundary conditions are 
used along the long (L) direction of the bar because they approximate the limit L — > oo 
better. We find that our results for the scaling function are different from those obtained 
for the film geometry and are also in reasonably good agreement with the most recent 
experimental results |13|] without any adjustable parameter. 

In order to describe the fluctuations of the order parameter in superfluid A He near the A 
point we employ the XY model which belongs to the same universality class as the Landau- 
Ginzburg free energy which is expressed in terms of the superfluid order parameter %p(r). 
Within the XY model the free-energy Tt can be written as 

Ti = -J^Si-sj, (3) 

where the summation is over all nearest neighbors, s = (cos 8, sin 8), and J > sets the 
energy scale. In this model Si are not real spins but they are pseudospin variables in 
which the angle 8 corresponds to the phase of the superfluid order parameter ip(r). The 
order parameter can be understood |L4|| as the average value of operator which creates an 



atom in the superfluid helium and it is defined in a volume whose linear extensions are 
much larger than the interparticle spacing and much smaller than the correlation length, a 
condition that is realized only very near the transition temperature. 

The pore geometry is represented by H x H x L lattices with L ^> H, namely, bar-like 
geometry. On the ends of the bar we always use periodic boundary conditions whereas we 
have carried out separate calculations for a) periodic and b) staggered boundary conditions 
along the sides (if-directions) of the bar. In the case of staggered boundary conditions the 
four sides of the bar are coupled to a staggered pseudospin configuration, i.e. the pseudospins 
at the boundaries alternate when we move along the lattice bonds (cf. also Ref. ||) so that 
the total sum of the spins at the sides of the bars is exactly zero. This boundary condition 
plays the role of the confining walls and corresponds to a boundary condition where the 
superfluid order parameter is exactly zero at the boundaries. In the simulation this can be 
easily visualized by corresponding every spin of the boundary layer to the same pseudospin 
but which interacts with all other spins of the layer next to the boundary via a coupling 
which alternates in sign. This spin is a dynamical variable which means that it can change 
direction during the cluster updating procedure and thus, the 0(2) invariance of the model 
is not broken at the boundaries. 

We computed the specific heat on H x H x L lattices, where H = 16, 20, 26, 32 and 
L = 5H. The specific heat c is obtained by 

c = ^ ((H 2 ) - (H) 2 ) , (4) 

where (3 = l/(ksT) and N is the number of pseudospins contributing to the specific heat. 
The multi-dimensional integrals in the expression for the averages in Eq. (f|) were computed 
by means of the Monte Carlo method using Wolff's 1-cluster algorithm fll5]| . To thermalize 
the system, we typically carried out of the order of 20, 000 steps involving single cluster 
flips during which no observable was calculated. After the end of that process we typically 
carried out of the order of 750, 000 single cluster flips where observables were calculated. 
The calculations were performed on a heterogeneous environment of computers including 
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IBM RS/6000 and DEC alpha AXP workstations and a Cray Y-MP. 
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FIG. 1. The specific heat as a function of t = T/T\ — 1 for 16 2 x L lattices with periodic 
boundary conditions (pbc) in all directions. In the three-dimensional XY model T\/J = 2.2017 
lh 



Here we check the finite-size scaling hypothesis for the specific heat with respect to the 
bar diameter H in the L — > oo limit. We do not need to take the actual L — > oo limit, 
because it turns out that lattices with the ratio L/H = 5 are a satisfactory approximation 
to that limit. We demonstrate this in Fig. [I| where the specific heat, calculated on lattices 
with ratios L/H = 5, 10, 15 and H = 16, is plotted near the peak. Within error bars the 
results for the three size lattices are the same. 

We are now able to compute the finite-size scaling function f\{tH x ' v ) defined by the 
expression 



c(t, H) = c(t , oo) + H^hitH 1 ^). 



(5) 



The function f\(x) is universal and v = 0.6705 as has been extracted from recent experi- 



ments fTF ]. The hyperscaling relation a = 2 — 3v yields a/u = —0.0172. At the reduced 
temperature t the correlation length £(t) = ^o = |t| _I/ becomes equal to the bar thickness H, 
i.e. t = (t$/H) 1 l v with $ = 0.498 ||. We have 
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c(* , oo) = c(0, oo) + 5+ to °, (6) 

where we use the bulk values c(0, oo) = 30, = —30 obtained by studying the finite-size 
scaling of the specific heat of cubes §| . 
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FIG. 2. Comparison of the finite-size scaling functions fiftH 1 ^) in Eq.([|) for films, bars and 
cubes with periodic boundary conditions (pbc) in all directions. 

In Fig. |^ we compare the scaling functions for cubes, films (taken from Ref. |5|) 

and bars in the presence of periodic boundary conditions. The more confining dimensions 
occur in the system, i.e. one, two, and three confining dimensions for the cube, the bar and 
the film, respectively, the more the function fi(x) is suppressed, only in the limit |x| — ► oo 
do the three functions agree. 

Here we compute the scaling function fi(x) again but using staggered boundary con- 
ditions in the if-directions of the lattice. These boundary conditions seem to mimic the 
physical boundary conditions met in the experiments rather well |||6|.|8| . Our results for the 
scaling function f\ are shown in Fig. [5] where fi is compared to that obtained with periodic 
boundary conditions. This figure demonstrates that in the presence of staggered boundary 
conditions the scaling function fi(x) is dramatically suppressed compared to the function 
fi (x) obtained for bar-like geometry with periodic boundary conditions in all directions. For 
large values of |x| both functions agree. Similar behavior was observed for a film geometry 
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FIG. 3. Comparison of the scaling functions fi(x) for bars with periodic and staggered bound- 
ary conditions on the sides. 



In order to relate our results to experimental results let us now compute the finite-size 
scaling function fi(x) obtained for bars with staggered boundary conditions in physical 
units. This leads to the same conversion formula used in Ref. fi(x)\ phys = A fi(x)\ lattice 
where the factor A is given by A = (Kti&b) / a 3 a) a ^ where V m is the molar volume of 4 He 
at saturated vapor pressure at T\. The unit of length a (i.e., the lattice spacing a in the XY 
model) was determined to be a = 2.95A |||§ and thus A = 15. 02 Joule/ (°Kmole). 

Coleman and Lipa ]13|] have recently measured the specific heat of superfluid helium 
confined in pore geometry. They have compared their results to the early data of Refs. 



nyi2fl. There is some range of disagreement near the critical temperature. The experimental 



results of Refs. (TT].[T2j have been obtained on much smaller diameter pores than those of Ref. 
fl"3| . In our recent work for the effect of the boundaries on superfluid films [K| we found that 
the boundary can create corrections to scaling which for films of thickness as large as those 
of the pore diameters of Refs. [0,0] can not be neglected. Since we want to stay away from 
such difficulties in this paper, we decided to use the results of the most recent work of Ref. 
jj~3| to compare with the results of our theoretical calculations. In Fig. |4] we compare the 
experimental results reported in Ref. 0] with our Monte Carlo results. Since the authors of 
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Ref . |L3| present only their specific heat data we deduced the scaling function fi (x) according 
to Eq.(|5|) using the bulk data of liquid 4 He reported in Ref. |TjJ. The agreement between the 
theoretical calculation and experiment is satisfactory. Note that the agreement between the 
function fx obtained from Monte Carlo data in the presence of periodic boundary conditions 
in the if-direction and the experimentally determined function fi would be far worse. For 
example, the peak of the function fi(x) for periodic boundary conditions in the if-direction 
in physical units is about 15Joule/(°Kmole) at x ~ —4, whereas for the experimental 
function fi(x) is approximately 6AJoule/(°Kmole)) at x ~ —23.1 (cf. Fig. §). The small 
difference between experimental and theoretical scaling functions obtained with staggered 
boundary conditions could be due to a number of reasons, including poor representation of 
the geometry or the boundary conditions in the real experimental system by the theoretical 
modeling. 
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FIG. 4. Our results for the specific heat scaling function fi(x) for pores (bars) with staggered 
boundary conditions on the substrate layers (open squares) are compared to the experimental 



results [13 1 (solid squares). For comparison the same function is also shown for the film geometry 
as obtained from our simulation (crosses, taken from Ref. ||) and from experiments (open circles, 
taken from Ref. H). H is expressed in units of A and fi(x) in Joule/(°Kmole). 



Fig. 0] demonstrates again the effect of the type of the confining geometry on the shape 



S 



of the universal scaling function fi(x), now in the presence of staggered boundary condi- 
tions which are a good approximation to the physical boundary conditions imposed by the 
confining walls. 

In conclusion, we have numerically computed the finite-size scaling function fi(x) for 
a bar geometry and we have compared our results to the scaling functions obtained earlier 
for cubes and films. For pores (represented by bars) we find good agreement between our 
Monte Carlo results in the presence of staggered boundary conditions and the most recent 
experimental results of Coleman and Lipa fl3||. More precise experimental data for the 
specific heat of liquid 4 He confined in pores near the A-critical temperature T\ are expected 
from future experiments [^0] which can be compared to our calculation. These experiments 
could employ different substrates to check the influence of the boundary conditions on the 
scaling function fi(x). Thus, this work together with our earlier work on other geometry 
|8|jl0[| shows that only scaling functions which belong to the same geometry and the same 
boundary conditions should be compared. 
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